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We present an efficient numerical code based on spectral methods to integrate the field equations of general 
Robinson-Trautmann spacetimes. The most natural basis functions for the spectral expansion of the metric 
functions are spherical harmonics. Using the values of appropriate combinations of the metric functions at 
the collocation points, we have managed to reduce expression swell when the number of spherical harmonics 
increases. Our numerical code runs with relatively little computational resources and the code tests have shown 
excellent accuracy and convergence. The code has been applied to situations of physical interest in the context 
of Robsinson-Trautmann geometries such as: perturbation of the exterior gravitational field of a spheroid of 
matter; perturbation of an initially boosted black hole; and the non-frontal collision of two Schwarzschild black 
holes. In dealing with these processes we have derived analytical lower and upper bounds on the velocity of the 
resulting black hole and the efficiency of the gravitational wave extraction, respectively. Numerical experiments 
were performed to determine the forms of the gravitational waves and the efficiency in each situation of interest. 



I. INTRODUCTION 

The use of numerical techniques to solve Einstein's 
field equations is a very promising strategy for dealing 
with problems of astrophysical interest such as gravi- 
tational collapse and formation of black holes, coales- 
cence of binary systems, the collision of black holes, etc. 
One could say that numerical relativity is the bridge be- 
tween relativity and astrophysics. A common feature in 
all these processes is the emission of gravitational waves 
- one of the most notable predictions of General Relativ- 
ity. Undoubtedly, detailed knowledge of how much sys- 
tem mass is converted into gravitational waves and the 
forms of the resulting gravitational waves are of funda- 
mental importance in the efforts to detect directly grav- 
itational waves with ground and space based observato- 



ries li. 

In this context we have two main objectives. The 
first is to present an efficient numerical code for solv- 
ing the field equations of general Robinson-Trautmann 
spacetimes, based on spectral methods ||2[]. The second 
is to exhibit the initial results for the wave forms and 
efficiency of the gravitational wave extraction associated 
with the following systems: the nonlinear perturbation of 
an oblate distribution of matter evolving to form a black 
hole, the recoil of a boosted black hole after interacting 
with a packet of gravitational waves, and the non-frontal 
collision of two Schwarzschild black holes. 

The Robinson-Trautman metric can be expressed 
as 131 
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where m is a null coordinate such that u = constant de- 



notes null hypersurfaces generated by the rays of the 
gravitational field, and that foliates the spacetime glob- 
ally; r is an affine parameter defined along the null 
geodesies determined by the vector d/dr and otq is a con- 
stant. The angular coordinates (9,(|)) span the spacelike 
surfaces u — constant, r = constant commonly known as 
the 2-surfaces. In the above expression an overdot indi- 
cates a derivative with respect to m, the functions A-, iden- 
tified as the Gaussian curvature of the 2-surfaces, and P 
depend on the coordinates m,9,(|). In the particular case 
of axial symmetry the metric functions do not depend on 
isf. Two important aspects of Robinson-Trautman space- 
times are worth mentioning: they are asymptotically flat 
and admit the presence of gravitational waves. 

Einstein's equations for the Robinson-Trautman (RT) 
spacetimes reduce to 



X = p2 + _^('sine^ 
sine V P 
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where the subscripts e,(|) denote derivatives with respect 
to the angles e and (|), respectively. The structure of the 
field equations is typical of the characteristic evolution 
scheme: Eq. (|2]l is the hypersurface equation that de- 
fines ^(M,e,(|)), and Eq. (O, known as the RT equation, 
governs the dynamics of the gravitational field. In other 
words, from the initial data F(Mo,e,(|)) the hypersurface 
equation determines A-(mo,9,'1'), and the RT equation al- 
lows us to evolve the initial data. 

Several works have been devoted to the evolution of 
the vacuum RT spacetimes focusing on the issue of the 
existence of solutions to the full nonlinear equation. The 
most general analysis on the existence and asymptotic 
behavior of the vacuum RT equation was given by Chr- 
usciel and Singleton ||4|], and independently by Frittelli 
and Moreschi Jst]. They proved that, for sufficiently 
smooth initial data, the spacetime exists globally for pos- 



itive retarded times, and converges asymptotically to the 
Schwarzschild metric. 

The combination of gravitational waves and the 
asymptotic convergence to the Schwarzschild solution 
allows RT spacetimes to be used to study the emission of 
gravitational waves in connection with the formation of a 
single black hole. The first work in this direction was due 
to Foster and Newman [6] who interpreted the linearized 
solution of the field equations as representing the emis- 
sion of gravitational radiation by a bounded source. RT 
spacetimes have also been used as test beds for numerical 
codes in the characteristic formulation of axisymmetric 
spacetimes 17|]. The full evolution of the field equations 
(|2]l and (O constitutes a valid and valuable framework for 
a detailed analysis of the wave forms and the efficiency 
of the gravitational wave extraction connected with the 
formation of black holes. Furthermore, a scenario in 
which the resulting black hole is moving with respect to 
an inertial observer due to, for example, the collision of 
two black holes, can be treated using the dynamics of 
RT spacetimes. Some of these aspects were addressed 
in Refs i8l 4l4ll in the context of axial symmetry after the 
development of a numerical code for integrating the field 
equations il5ll . 

The paper is organized as follows. In Section 2 we 
present the details of our numerical scheme. In Section 3 
we introduce the initial data employed for the code tests. 
Tests confirm the excellent accuracy and convergence of 
the numerical code. Section 4 is devoted to applying the 
code to situations of physical interest, with the gravita- 
tional wave forms and the efficiency of the gravitational 
wave extraction being evaluated after evolving the space- 
time. In addition, by analyzing the conservation of the 
four-momentum, we were able to determine analytically 
a lower limit of the velocity of the black hole formed as 
well as an upper bound for the efficiency of the gravita- 
tional wave extraction. As a final and brief application, 
we consider initial data describing the collision of two 



Schwarzchild black holes in three situations: orthogonal, 
oblique and head-on collisions. Finally, in Section 5 we 
present a summary of our results and conclusions. 

II. THE NUMERICAL METHOD 



The modes (ai^j^bi^j) are determined such that the resid- 
ual vanish in the mean jlSll . This means that the inner 
products (projections) of the residuals with respect to 
the test functions are made equal to zero. For the G-NI 
method the test and basis functions are the same, then 



According to the Galerkin method with numerical in- 
tegration (G-NI) il6ll . the metric functions f (m,9,(|)) and 
A-(m,9,(|)) are approximated by a series expansions given 
by, 

N l=k 

Pn{u,BA) = Y. L ''kM^kiiQA) (4) 

k=0 l=-k 
N l=k 

XNiu,eA) = Y. L ^kiiu)Yki{eA), (5) 

k=Q l=-k 

where yjtj (9 , (|)) are the spherical harmonics chosen as the 
basis functions, ^^/(m) and /?,(./(«) are the unknown modes 
and N is the truncation order that dictates the number of 
independent modes. These modes must satisfy the sym- 
metry relations a^i = (— l)^'fl,t/ and b^i ~ {—l)^'bki 
that result in a total number of 2(A^+ 1)^ independent 
modes. Although the number of independent modes in- 
creases considerably with the truncation order, the use of 
spherical harmonics provides exponential convergence 
for smooth functions on the sphere, and also fixes the 
so called pole problem 111711 . 

The substitution of the approximations Eqs. dH and © 
into the field equations (|2|i and Q results in the residuals 
associated to these equations, 

Resx{u,d,^)=XM-PNY,{k^ + k-l)auYu + 
k,l 
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where Res(u,9,(|)) stands for Resx(M,9,(|)) or 
Res/>(M,9,(|)), with m varying from to A^, and 
n = — m, ..,m. 

As dictated by the G-NI method, these integrals were 
evaluated using quadrature formulas which provide fast 
and accurate results, since these formulas usually con- 
verge exponentially 1171 Il9ll . We introduce the variable 
X = cos 9 and define the quadrature collocation points 
{xk,(^j)hy 
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^j = f J = 0,l,..,M 
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xii = —I, zeros of !P~ (x) , 1 , 
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where iP^(x) is the Legendre polynomial of order A^, and 
we have set M = 2N and A^ = 3N/2 that leads to a more 
accurate integration. The integration with respect to (|) 
the trapezoidal rule was used that is the indicated quadra- 
ture rule for periodic functions 1 17], while the Legendre- 
Gauss-Lobatto quadrature formula was used for the inte- 
gration with respect to x. These integrals can be written 
schematically as 



M N 



/ / F{x,^)dxd^^yyvjWkF{xk,^j) (10) 

where Wk,Vj are the weights 11911 . 

In order to allow the use of high truncation orders, the 
projections of the residual equations (|6]l and Q are ex- 
pressed as, 



bkj = Skj aimPl,m{P.Q)lmi 
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where P = Pn and the indices Im indicate that the expres- 
sion is evaluated at the collocation points x = x/,(|) = (|)„,. 
These values are related to the modes Okj but stored 
in separated files. Also, A-/„, represent the values of 
X{u,Q,(^) ='E.kj k{^+^)bki{u)Yki{Q,<^) at the collocation 
points. Therefore, the field equations (O and (|3]l are re- 
duced, respectively, to a set of algebraic equations con- 
necting the modes bjk,aii„, and a set of ordinary dif- 
ferential equations for the modes a/,,,. The evolution 
scheme proceeds as follows: (i) from the initial data 
Po{x,(^) — P{uo,x,(^) the initial modes ajk{uo) are fixed 
from Eq. (|4|i; (ii) the values of P,Pq and P(j,/sin9 at the 
collocation points are evaluated; (iii) the initial modes 
bjk{uo) are determined from the set of equations ( fTTl i: 
(iv) the values of A- at the collocation points are evalu- 
ated; (v) the ordinary differential equations (fT2] | fix Ojk at 
u = UQ, and the modes Ujk can be determined at the next 
null surface uq + 5m. By repeating the whole process the 
evolution of the modes is determined, and consequently 
the evolution of the functions P{u,x,<^) and X{u,x,(^). 

III. CODE TESTS 

The integration of the dynamical system (fT2] ) re- 
quires (A^+1)^ initial conditions akj{uo), determined 
by the initial data Po(-^,^) — P{uq,x,(^). We have de- 
termined physically relevant initial data for the axisym- 
metric case representing the exterior gravitational field 
of homogeneous and non-homogeneous spheroids lll2ll . 
the perturbed boosted black hole by a gravitational wave 
packet Ollll and also the collision of two Schwarzschild 
black holes OlOll . We consider here the generahzation 
of these initial data to the non-axisymmetric case, but 



details of the derivation will be omitted since they are 
straightforward. 

We start with simple initial data which have no direct 
physical motivation. 



'^,i'^(^,^) 



l-e2(i_jc2)sin2(|) 
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where < e < 1 plays the role of the eccentricity of the 
above spheroid-type function. The second set of initial 
data represents the exterior gravitational field of a per- 
turbed oblate spheroid lll2ll . 
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where the first two terms on the right describe a homoge- 
neous oblate spheroid in which Bq is an arbitrary param- 
eter, a = arctan(l/Q, p = (1 +3Z,^)aTctSLn{l/Q, with ^ 
being a free parameter associated with the oblateness of 
the spheroid 1121 . The function /(x,(|)) can be interpreted 
as describing a perturbation inside the matter distribution 
of the spheroid. The third set of initial data is expressed 
as 
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(x,(j)) 



1 




•\/Po(coshy±xsinhy) 

where Pq and y are arbitrary parameters, and g{x,<^) is 
a regular function. This function generalizes the initial 
data of Ref. 11 ill that describes a boosted black hole with 
respect to an asymptotic observer which is perturbed by a 
non-spherical mass distribution. The boosted black hole 
in the z— direction is characterized by the exact stationary 
solution 12011 of the field equations Q and ([3]l, 



P = Po(coshY±jt:sinhY), 



(16) 



in which X = Pq. The black hole is moving with velocity 



tanhy with respect to an asymptotic observer ll20ll . 



This solution has been used to study the deceleration of 
a moving black hole due to an external perturbation and 
also the head-on collision of two black holes with the 



same velocity lIlOll . 

Due to the lack of exact non-stationary solutions of 
the field equations, the best way of testing the accuracy 
and convergence of the code is to verify if the conserved 
quantity 02111 



1 
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is maintained constant by the numerical solution within 
an acceptable error, and if the error diminishes when 
the truncation order A^ increases. Since the initial data 
Pq{x,(^) determine the exact value, / = /exact, a useful 
measure of the numerical error with respect to the exact 
value of / is 



5/ 



(18) 



T ' 

^ exact 

where /numer is the numerical value of/. 

In our numerical experiments we have determined the 
maximum deviations after evolving each family of ini- 
tial data for several values of A^ as indicated in Fig. 1. 
In the first family of initial data ( fT3T l we set e = 0.6 
(note that a very small value for e results in a weak per- 
turbation of the stationary solution, which corresponds 
to e = 0). For the second initial data (fT4l l we have 
set Bo ~ 0.6,^ = 0.1 and used a Gaussian-like pertur- 
bation /(jc,(^) = Aosin2(j)(l-x2)2e-('^-''-^)', with Aq = 
0.05. The third family of initial data ( fTSI l is used with 
Po = 1.0,y = 0.5, while the perturbation is given by 
g{x,(if) = Ao(l +sin(j))(l -x2)2e-(-'^-oi)', with Aq = 0.1. 
In Fig. 1 the log-linear plots of the maximum relative 
error versus the truncation order A^ are presented for all 
three cases. Increasing A' produces an exponential de- 
crease of 5/ until the saturation due to the round-off error 
is achieved. 
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FIG. 1. Plots of the maximum relative errors 8/ after evolv- 
ing the initial data | |13I ). ( I14t and | |15I ) (from top to bottom). 
In all cases the exponential decay of the maximum error was 
observed as a result of the fast convergence of the code. The 
saturation due to the round-off error is shown for the two first 
cases and is achieved for different values of A'. 

Another code test is to evaluate the L2 norm associated 
with the residual evolution equation (|7]l given by. 



Liiu) 



2% 



Res^(M,x,(|))(i(|)(ijt: 



1/2 



(19) 
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We have followed the evolution of L2 for each set of ini- 
tial data and several truncation orders and selected the 
corresponding maximum values. The results are shown 
in Fig. 2 and show spectral convergence. 
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FIG. 2. Plots of the maximum values of the L2 norm (see Eq. 
J19H with respect to the truncation orders A' for the initial data 
il3t . l ll4t and l ll5t (from top to bottom). Again, spectral con- 
vergence is obtained. 



The third code test is to reproduce correctly the 
Schwarzschild solution as the asymptotic configuration 
of RT spacetimes. In the axisymmetric case the asymp- 
totic configuration can be either a black hole at rest or 
moving with constant velocity along the symmetry axis 
with respect to an inertial observer The boost is a con- 



sequence of the total energy-momentum conservation 
when the flux associated with the gravitational waves is 
nonzero. The net momentum of the gravitational waves 
is transferred to the source producing a boosted black 
hole along the symmetry axis. In the general case we 
expect that the resulting black hole will be moving in an 
arbitrary direction. 

The structure of the initial data plays a crucial role in 
determining whether or not the black hole acquires mo- 
mentum. Let us evolve the field equations taking the first 
set of initial data (Eq. ( fOl l) with several values of the 
parameter e. In all cases the aoo mode tends to a constant 
value while all other modes approach zero. The final con- 
figuration is Pflnai = floo = constant, which is iden tified a s 
the Schwarzschild solution with mass-energy ||8, 
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In fact, the absence of any boost (y = 0) is a consequence 
of the symmetric pattern of gravitational wave radiation 
produced by the initial configuration (see Fig. 3). 

We turn now to the less symmetric second and third 
sets of initial data. After several numerical experiments 
we have obtained the following asymptotic expression 
for the metric function P, 



Pfinai (jc, (j)) = ao -I- aix + a2 vl-^cos (|) -I- 

aaVl-x^sinc]) (21) 

where a^, y = 0, 1 , .., 3 are determined numerically. This 
expression generalizes the stationary solution ( fTSI l, rep- 
resenting a boosted black hole in an arbitrary direction 
with respect to an asymptotic observer as shown by Cor- 
nish 02311 . The parameters Uj are not independent but 
satisfy the constraint ttg — L/=i ^C/ = ^0 = constant im- 
posed by Eq. (|2]l. This allows the following alterna- 
tive choice for these parameters: ao = Po coshy, ai = 
Poflsinhy, a2 = Pobsinhj, aj, = Pocsinhy, which yields 



sinhyvl — x2(Z7Cos(|) + csin(|))], (22) 



and the relation, 
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(23) 



that must be satisfied by the numerical solution asymp- 
totically. Notice that b = c = implies a ~ ±1 recover- 
ing the axially boosted black hole (fTSI l. where the signs 
± indicate opposing directions of movement along the z- 
axis. In the general case, the black hole is moving with 
respect to an inertial observer in an arbitrary direction 
determined by a,b , c, which can be identified as the di- 
rection cosines lll3ll . In this case the total mass-energy of 
the black hole is given by. 



Ebh = -r- 
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where mrest = nioPQ is the rest mass of the boosted 
black hole. The velocity of the hole is thus given by 
tanhy as in the axisymmetric case. We have checked 
the relation (|23]) with the asymptotic numerical solution 
obtained after evolving the third family of initial data. 
The results are shown in Table I with the deviation of 
(a^ + b^ + c^)numer from the exact value versus the trun- 
cation order A^. Accordingly, for the smallest truncation 
order, N = 6, the deviation is about one part in 10^^ con- 
firming the excellent accuracy of the code. 

IV. WAVE FORMS AND EFFICIENCY OF 

GRAVITATIONAL WAVE EXTRACTION. 

NON-FRONTAL COLLISION OF TWO BLACK HOLES 



One of the most important aspects of RT spacetimes 
is the presence of g ravitational waves |3]. As discussed 
in Refs. ialloHull, the characterization of gravitational 



TABLE I. The deviation | (a^ + b^ + c^)„ 
increasing truncation orders. 



1 1 evaluated for 



waves in RT spacetimes is based on the Peeling the- 
orem [24], for which the Weyl tensor is expressed as 
Cabcd ^ NABCo/f at large r characterizing the wave 
zone. The relevant components of the Weyl tensor 
with respect to a semi-null tetrad basis tSJ are C0303 = 
-C0202 = -£>(«, -J^, (]))/'•+ o(l/r2), where 



1 Pm 



(25) 



(l-x2) P 

This function gives the time and angular dependence of 
the gravitational wave amplitude within the wave zone. 

The structure and evolution of the angular pattern of 
gravitational waves is directly connected to the initial 
data. In the next numerical experiments we show the 
evolution of these patterns for different initial data by 
displaying a sequence of three dimensional polar plots of 
Z)(m,9,(|)). In Fig. 3 the first and second sequences corre- 
spond to the initial data (fT3] i and dTsl l, respectively, where 
D(m,0,(|)) is evaluated at several times. The symmetry in 
the angular distribution of D in the first sequence implies 
that no momentum from the gravitational waves is trans- 
ferred to the source. A distinct situation occurs for the 
second sequence in which a non-symmetric pattern of D 
indicates the existence of a net flux of momentum carried 
away by gravitational waves. In this case, due to the con- 
servation of momentum-energy, the end state is a boosted 
black hole described by Eq. ( l22b . Another example of a 
non-symmetric angular pattern is shown by a sequence 



of two-dimensional plots in Fig. 5. Note the presence of resents an oblique collision of two Schwarzschild black 

lobes pointing in the directions of maximum emission of holes (cf. Eq. ( |33] |) whose end state is, in general, a 

gravitational radiation. In this case the initial data rep- boosted black hole. 

I 




♦ ♦ ♦ 







FIG. 3. Sequence of three dimensional polar plots of D(u,Q,^) corresponding to the first set of initial data evaluated at u = 
0,0.1,0.4,0.7, and to the third set of initial data evaluated at u = 0.002,0.004,0.006,0.008. The symmetric pattern in the first 
sequence suggests that no momentum is transferred to the remnant, while in the second sequence the lack of symmetry of the 
angular pattern indicates that the resulting black hole might be moving with respect to an asymptotic observer. 



The efficiency of gravitational wave emission mea- 
sures the amount of mass-energy extracted during the 
evolution of RT spacetime with the eventual formation 
of a black hole. It is expressed as 



A= = 1 

El 



miestcoshy/ 
El ' 



(26) 



where Ei is the mass-energy of the initial system, Ebh = 
TOiestCOshy/ is the energy of the remnant black hole and 
Y/ is the final boost parameter. We can obtain analytical 
information about the efficiency of any process described 
by general Robinson-Trautmann spacetimes by analyz- 
ing the conservation of global four-momentum. Follow- 
ing Eardley ll25ll we have. 



El — Ebh + Eqw 



Pi = Pbh + Pgv/, 



(27) 



(28) 



where iicw and pcw are the amount of energy and mo- 
mentum carried by gravitational waves; pi and pbh are 
the initial and black hole momentum magnitudes. These 
quantities are evaluated from the components of the four- 
momentum £P" associated with RT spacetimes 1 1 3i J 
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where /" = (1, \/T~~?cos(|),\/T~^x2sin(|),x). Since 
the remnant is a black hole, it follows that pbh = 
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FIG. 4. Illustration of the three situations representing (from 
top to bottom) orthogonal, oblique and head-on collisions of 
two Schwarzschild black holes. 

mrestsinhy/; also iscw > PGW as a consequence of the 
radiated four-momentum being future pointing and time- 
like. This inequality is conveniently re-expressed using 
Eqs. ^ and ^ as, 




FIG. 5. Illustration of the sequence of polar plots of D pro- 
jected on the planes x = (m = 0.02,0.10,0.26) and x = 0.3 
(i( = 0.02,0.10,0.18) for the oblique collision. The angular 
pattern is dominated by lobes which indicate the directions of 
maximum magnitude of the gravitational wave emission. Note 
that the two dominant lobes open as the black hole is deceler- 
ated, forming a bremsstrahlung-like pattern. 
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coshy/- — sinhy/^ < = p > 0, 

'Wrest 



(30) 



since to E\ > pi. Now, from the initial data P{uo,x,(^) 
Ei,pi can be evaluated via Eq. ( |29] i; also njrest is de- 
termined from the conserved quantity / [j. As a conse- 



quence, P becomes determined after specifying the initial 
data. A direct manipulation of the above expression re- 
sults in the determination of the lower bound of the boost 
parameter, or Yf > Ymin — arccosh[(p^ + 1)/2P]. In this 
way, the lower bound of the remnant velocity Vf = tanh Yf , 
is given by, 



By calculating the conserved quantity / from Eq. I17t with P given 
by Eq. (22) we found / = Pq ■ However, since / has its value fixed 
by the initial data, it follows that /"o = ^ - which determines the rest 



mass of the resulting black hole as m„ 



- mo/Pfl 
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(31) 



From the lower bound of the final velocity it is possible to 



estimate the upper bound of the efficiency for any initial 
data. 



^max — t 
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As a brief appHcation of these estimates, we study 
the efficiency of the non-frontal collision of two 
Schwarzschild black holes initially with the same speed 
V = tanhy described by the following initial data, 

P(M0,e,^)=(^=4== + ^=4==| (33) 



Vi\(m vpW, 



where 



P± == pf^[coshy 



fl±sinhycos9 



sinhysin0(/7±cos(|) + c±sin(|))] with a\+b\+c\ = 1. 
Note that the particular case b± — c± — Q, a± — ±1 
reduces to the initial data for two initially boosted 
Scharwzschild black holes with opposite velocities 
tanhy along the symmetry axis llOj]. In the general 
case the black holes have initially arbitrary directions 
determined by the parameters a,b,c. By conveniently 
setting these parameters we can study three situa- 
tions illustrated by Fig. 4, that is, two non-frontal 
collisions and one frontal collision. In all situations 



10 

the initial motion of the first black hole is along the 
z-axis, with «+ = 1,/?+ = c+ = 0, whereas the sec- 
ond black hole has in turn (i) a_ = /?_ = 0,c_ = 1 
corresponding to a direction orthogonal to the z-axis, 
(ii) fl_ = -V2/2,b^ = 0,c_ = -\/2/2, and (iii) 
fl_ = — !,/?_ = c_ =0 for a head-on collision. In all 
numerical experiments the boost parameter was fixed 
as y = 0.9 which corresponds to the initial relativistic 
velocity v « 0.716, P^ = 1 and Pi is a free parameter 
that defines the ratio r| between the initial rest masses of 
the two black holes. 



'rest ^ IpiO) 



(34) 

In table II we summarize the results for several values of 
r|. Note that the frontal collision is the most efficient in 
terms of carrying away energy by gravitational waves. In 
line with numerical simulations based on more realistic 
scenarios, the efficiency of the extraction of energy via 
gravitational waves never exceeds a few percent of the 
initial energy. The results are also consistent with analyt- 
ical estimates for both the lower bound of the final boost 
parameter and the maximum efficiency, as given by Eqs. 
dJTT i and ( l32b respectively. The values of the maximum 
boost parameter are displayed in Table II for each case. 





Orthogonal 


Oblique 


Head-on 


n 


^'min 


^'f 


A(%) 


^'min 


^•f 


A(%) 


^min 


^'f 


A(%) 


0.001 


0.597 


0.601 


0.608 


0.525 


0.537 


1.23 


0.498 


0.514 


1.473 


0.008 


0.557 


0.564 


0.911 


0.437 


0.456 


1.97 


0.383 


0.411 


2.428 


0.125 


0.509 


0.520 


1.307 


0.306 


0.327 


3.019 


0.160 


0.202 


3.797 


0.296 


0.501 


0.511 


1.384 


0.269 


0.294 


3.233 


0.077 


0.121 


4.084 



TABLE II. Efficiency and the final boost of the remnant black hole for several values of Tj. 



V. DISCUSSION 

In this paper we have studied numerically the dynam- 
ics of general Robinson-Trautman geometries. We have 



I 

constructed an efficient and accurate numerical code 
based on the Galerkin method that generalizes the code 
we have implemented for the axisymmetric case 111511 . 
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Spherical harmonics were used as the basis for the spec- 
tral expansions of the metric functions P(m,9,(|)) and 
A,(m,9,(|)) (cf. Eqs. (|4| and Q, respectively). How- 
ever, as can be seen from these expansions, an increase 
in the truncation orders A^ and M increases consider- 
ably the number of modes akj^bkj- The residual equa- 
tions associated with the field equations were projected 
using the quadrature formulas Eq. ^. The field equa- 
tions (|2]i and ^ were reduced to a set of ordinary dif- 
ferential equations for the (A^+ 1)^ modes a^j together 
with equations for the modes bkj, the values of A,, P 
and their derivatives (see Eqs. (fTTI) and ( fT2] l) at the col- 
location points. This procedure allowed us to reach, 
without much computational effort, a truncation order 
of A^max = 14, corresponding to 225 modes ai^j and an 
equal number of independent modes bkj- The code tests 
consisted of checking the conserved quantity I, monitor- 
ing the L2-norm associated with the residual RT equa- 
tion, and reproducing the correct asymptotic solution de- 
scribed by Eq. ( 1211 1 of a black hole moving with constant 
velocity in a arbitrary direction with respect to an asymp- 
totic observer 12311 . These tests confirmed the excellent 
accuracy and exponential convergence of the code. 

Another important result is the generalization of phys- 
ically relevant initial data deduced for the axisymmet- 



ric case lllOl4l2ll . Though we have not presented the de- 
tails of these derivations here, we have established initial 
data representing the exterior gravitational field of a per- 
turbed oblate (prolate) spheroid of matter, for which the 
perturbation can be a disruption inside the source. We 
have also generaUzed the initial data representing a per- 
turbed axisymmetric boosted black hole. In both cases 
the asymptotic solution is a boosted black hole in an ar- 



bitrary direction determined by the parameters a^b,c (cf. 
Eq. ^). 

We have presented two applications of the code. The 
first is a visualization of the angular pattern of the gravi- 
tational radiation defined in the wave zone (cf. Eq. (l25Tl). 
whose structure depends on the initial data. The second 
concentrates on the efficiency of the conversion of the 
initial mass to gravitational waves. In this case we have 
set up the equations governing the conservation of global 
four-momentum taking into account that the remnant is a 
boosted black hole moving in an arbitrary direction with 
respect to an asymptotic observer. Analytical estimates 
of the minimum or lower bound of the boost parameter 
and the maximum efficiency were derived from the initial 
data. 

Finally we point out that a thorough study of the non- 
frontal collision of two Schwarzschild black holes in the 
realm of RT spacetimes is a subject for further investi- 
gation. The initial data describing a non-frontal collision 
of two Schwarzschild black holes with the same speed, 
given by Eq. ( l33T l, generalizes the initial data established 
for the head-on collision presented in Ref. 1101 . The 
numerical integration of the field equations (|2]| and (O 
provides a useful theoretical probe for studying the con- 
sequences of the non-frontal collision of two black holes 
02611 . In our preliminary analysis we have confirmed that 
a head-on collision is the most efficient source for ex- 
tracting gravitational radiation, although obviously be- 
ing a rarer astrophysical event than a non-frontal colli- 
sion. The evolution of the apparent horizon would also 
be worth investigating. 

The authors acknowledge the financial support of the 
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